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Abstract. A very rich phase diagram has recently been found in CeCu2Si2 from high pressure experiments 
where, in particular, a transition between an intermediate valence configuration and an integral valent 
heavy fermion state has been observed. We show that such a valence transition can be understood in the 
framework of the periodic Anderson model. In particular, our results show a breakdown of a mixed-valence 
state which is accompanied by a drastic change in the / occupation in agreement with experiment. This 
valence transition can possibly be interpreted as a collapse of the large Fermi surface of the heavy fermion 
state which incorporates not only the conduction electrons but also the localized / electrons. The theoretical 
approach used in this paper is based on the novel projector-based renormalization method (PRM). With 
respect to the periodic Anderson model, the method was before only employed in combination with the 
basic approximations of the well-known slave-boson mean-field theory. In this paper, the PRM treatment 
is performed in a more sophisticated manner where both mixed as well as integral valent solutions have 
been obtained. Furthermore, we argue that the presented PRM approach might be a promising starting 
point to study the competing interactions in CeCu2Si2 and related compounds. 

PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) - 71.27.+a Strongly correlated electron 
systems; heavy fermions - 75. 30. Mb Valence fluctuation, Kondo lattice, and heavy-fermion phenomena 



1 Introduction 

Since the discovery pQ of a superconducting state formed 
by heavy quasi-particles in CeCu2Si2 this and related com- 
pounds have attracted a lot of scientific interest. Despite 
its long history only recently a whole variety of new physi- 
cal phases has been observed which was possible by the in- 
triguing development of experimental techniques. By sub- 
stituting Si by Ge in the parent compound CeCu2Si2 a 
continuous change from a heavy fermion (HF) supercon- 
ducting phase to an antiferromagentic state was observed 
|2] . An even more complex phase diagram has been found 
in pure CeCu2Si2 by applying high pressure There 
two superconducting phases with different pairing mech- 
anisms have been found besides an antiferromagnetic and 
a HF phase. Furthermore, a transition between interme- 
diate and integral valence states has been observed. (For 
a recent review on superconductivity in Ce based HF ma- 
terials see Ref. 0.) 

From the theoretical point of view the periodic Ander- 
son model (PAM) is considered to be the basic microscopic 
model for the investigation of HF systems |Bj • The PAM 
describes the interaction between localized, strongly cor- 
related / and itinerant conduction electrons. In the limit 
of infinitely large Coulomb repulsion on / sites the PAM 



can be written as 



H — Ho + Hi, 



Ua - £ f E hmfim + ^2 Ek C km C km , 
i 7 rn k.m 

^ = ^ E ^(/Lc kme ikR < + h. c .) 



(i) 



Here, £f and £k> both measured from the chemical poten- 
tial, are the excitation energies of localized / and itiner- 
ant conduction electrons. As a simplification, often both 
types of electrons are assumed to have the same angular 
momentum index m with Uf values, m = \...Uf. The in- 
finitely large local Coulomb repulsion is taken into account 
by Hubbard operators 

flm = fim IT — fimfi™) 

which enable either empty or singly occupied / sites. 

Due to the complexity of the PAM, most theoretical 
studies only focus on certain aspects of the rich phase dia- 
grams of rare earth materials. Slave-boson mean-field (SB) 
methods, large-N expansions, and the dynamical mean- 
field theory [7] have been applied to discuss the interplay 



2 



A. Hiibsch and K. W. Becker: Valence transition in the periodic Anderson model 



between RKKY and Kondo interactions. Thereby, a tran- 
sition between an antiferromagnetic phase and a param- 
agnetic state was discussed. On the other hand, to de- 
scribe the valence transition and HF superconductivity in 
CeCu2Si2 an extended PAM was studied. This model in- 
cludes an additional Coulomb interaction between / and 
conduction electrons and was discussed within a slave- 
boson fluctuation approximation 

In this paper we apply a novel projector-based renor- 
malization method (PRM) J§] to the PAM with the aim to 
address the question whether a valence transition, as ex- 
perimentally observed in CeCu2Si2 [30, can occur in the 
plain model. For that purpose we extend in this paper our 
previous work on the PAM JHj, which was restricted to 
the HF phase. The PRM provides a natural way to discuss 
the interplay of competing interactions which naturally 
emerge from the renormalization treatment of the PAM. 
Therefore, we believe that the PRM represents a suited 
approach for a deeper understanding of the rich phase di- 
agram of CeCu2Si2 or of related compounds. However, in 
this paper we concentrate on the valence transition, never- 
theless, we are able to sketch how superconducting phases 
and RKKY interactions could also be included in our ap- 
proach. 

This paper is organized as follows. In the next section 
we briefly describe the novel PRM approach [Hj that is 
applied to the PAM in Sec. 01 Here, the Hubbard opera- 
tors, introduced to take into account the infinitely large 
Coulomb repulsion on / sites, cause the main problems of 
any theoretical approach. It will turn out that the well- 
known SB theory |lllll2j as well as our recent analytical 
approach based on the PRM do not sufficiently pre- 
vent from unphysical states with doubly occupied / sites. 
In contrast, the modified PRM treatment of Sec. |3| strictly 
suppresses doubly occupied / sites by taking into account 
electronic correlations by means of the Hubbard opera- 
tors. Results are presented in Sec. 0] where mixed valent 
as well as integral valent states are found, and a valence 
transition is observed. Furthermore, wc compare our re- 
sults with the solutions of the SB theory and our PRM 
approach of Ref. Finally, we summarize in Sec. El 



2 Methodology 

The PRM approach [5] starts from a decomposition of a 
given many-particle Hamiltonian, H — Ho + Hi, where 
the perturbation Hi should not contain any terms that 
commute with the unperturbed part Ho- Thus, Hi repre- 
sents transitions between eigenstates of Ho with different 
eigenenergies. In the following, we assume that the eigen- 
value problem of Ho is solved, 

H \nW) = £#V°>). 



A crucial idea of the PRM is the definition of projection 
operators by 

P X A = l» (0) > (™ (0) I <™ (0) l-4|™ (0) > (2) 
x0(A-|4°)-^°)|). 

Note that Pa and Qa = 1 — Pa are super-operators acting 
on ordinary operators A of the unitary space. Pa projects 
on those parts of A which are formed by transition opera- 
tors |n(°))(m(°)| with energy differences \En^ — Em \ less 
than a given cutoff A. (A is smaller than the cutoff A of 
the original model.) On the other hand, Qa projects on 
the high-energy transitions of an operator. Note, in par- 
ticular, that in Eq. ijSJ) neither |n^ ^) nor \m^) have to be 
low-energy eigenstates of Hq. 

Next, an effective Hamiltonian H\ is derived from the 
original Hamiltonian H by an unitary transformation 

Hx = e x >- H e~ Xx , (3) 

where the anti-Hermitian generator of the transformation, 
X\ = — X^, shall be chosen in such a way that only tran- 
sition operators (between eigenstates of Hq) with transi- 
tion energies less than the given cutoff A contribute to H\ . 
Thus, the condition 

QxHx = (4) 

must be fulfilled and will be used below to determine X\ . 
Note that it is straightforward to evaluate Eqs. J3J and (0} 
in perturbation theory ||. However, using an appropriate 
ansatz for the generator X\, the effective Hamiltonian H\ 
can also be calculated in non-perturbative manner. 

A renormalization scheme can be derived if the elimi- 
nation procedure for the interaction Hi is not performed 
in one step but rather a sequence of unitary transforma- 
tions of the form 

W(A-zu) =e x ^Hxe~ x ^ (5) 

is applied to the original Hamiltonian H. Thus, transi- 
tions between eigenstates of Hq caused by the interaction 
Hi are eliminated in steps where the respective transition 
energies are used as renormalization parameter A. Fur- 
thermore, 

Q(\-A\)H(X-A\) = (6) 

is used to specify the generator X\^a\ of the unitary trans- 
formation. Note that Eqs. JHJ and © describe a renormal- 
ization step that decreases the cutoff of the Hamiltonian 
from A to (A — AX), as one can see from a comparison with 
Eqs. J2J) and (@J. Therefore, difference equations for the A 
dependence of the Hamiltonian can be derived from (JSJ 
and and we call the resulting equations for the pa- 
rameters of the Hamiltonian renormalization equations. 
Note, that the solutions of these renormalization equa- 
tions strongly depend on the parameters of the original 
Hamiltonian H, and that the limit A — > provides the 
desired effective Hamiltonian without any interactions. 
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3 Renormalization of the PAM 

In the following, we want to apply the framework of the 
PRM as discussed above to the PAM. It is well known 
that much of the physics of the PAM JIJ can be un- 
derstood in terms of an uncorrelated model, that is for 
vanishing Coulomb repulsion on / sites where the Hub- 
bard operators f} m are replaced by usual fermionic opera- 
tors j\ m . This model can be solved exactly. However, the 
parameters have to be renormalized appropriately. Var- 
ious theoretical methods have been developed to gener- 
ate renormalized Hamiltonians. Most popular is the one 
derived from slave-boson mean-field (SB) theory |11I12| . 
Note however that only HF type solutions can be obtained 
in this way. In particular, the SB solution breaks down if 
the original / level is located too far below the Fermi level 
or, equivalently, if the hybridization strength between 
/ and conduction electrons becomes too weak. 

In the Hamiltonian the Hubbard operators f\ take 
care of the infinitely large local Coulomb repulsion on / 
sites so that multiple occupied / sites are strictly forbid- 
den. Consequently, any effective model has to satisfy this 
requirement as well. However, SB mean-field theory as 
well as our recent PRM treatment of the PAM PHj map 
the Hamiltonian of Eq. QJ onto an effectively free system 
consisting of two non-interacting fermionic quasi-particles, 

Weff = Y W k 4m C km + Y W k ftmhm + ^shift, (7) 
k,m k,m 

with renormalized parameters w£, ujI, -E s hift- It is impor- 
tant to notice that due to construction the effective Hamil- 
tonian Tl c g of Eq.Q does not prevent from multiple oc- 
cupation of / sites. This follows from the occurrence of 
the fermionic operators /k m and f^ m in Eq.(0) instead 

of the Hubbard operators /k m and f^ m - However, an ap- 
proximation that involves a replacement of the Hubbard 
operators /km and f^ m by usual fermionic operators /k TO 
and fJ might lead to useful results as long as only very 
few / type states are below the Fermi level. Thus, only 
HF-like solutions with a renormalized / level above the 
Fermi level can be obtained based on effective Hamiltoni- 
ans of type 0, and SB mean- field theory as well as our 
recent PRM treatment of the PAM JOj can not describe 
integral valent states. 



3.1 Renormalization ansatz 

In this paper we want to describe the transition of the 
PAM between mixed valent and integral valent states. 
Thus, a theoretical treatment is needed that reliably pre- 
vents from unphysical multiple occupation of / sites. For 
that purpose, we again apply the framework of the PRM 
to the PAM, but, in contrast to our recent work JOj, we 
now keep the Hubbard operators during the whole renor- 
malization procedure. Thus, the renormalization ansatz 



reads 

7~(\ = 7~Lo,\ + 7~Li,\ (8) 

H ,A = fif.X Y /km/km + Y ^k,A (/L/kra) NL 
k , in k,m 

+ Y £k > A C km C km + E X 
k , in 

H t ,x = Pa^i = £V k Pa (/L c k- + h.c.) 

k,m 

after all excitations with energies larger than the cut-off A 
have been eliminated. Due to the renormalization process 
all parameters depend on A, and an additional energy shift 
E\ and a hopping between different / sites, 

f ft f \ - J_ V P I e ik ( R '- R ^ 

i,j(¥*) 

have been generated. Here, we have introduced Fourier 
transformed / operators, 

ft = _L V P e lk R * 

J km v/jV ' lm 

The initial parameter values of the original model (at cut- 
off A = A) are 

Pf,A = £f, ^k,A = 0, e k ,yi = £k, Ea = 0. (9) 

To perform the PRM scheme we also need the commu- 
tator of the unperturbed Hamiltonian with the hybridiza- 
tion. For convenience, we introduce the unperturbed Liou- 
ville operator Lo,a which is defined by IjqxA — [Ho,X>-A] 
for any operator variable A, and to simplify the calcula- 
tions, the one-particle operators and c£ are consid- 
ered as approximate eigenoperators of Lo,a, 

Lo,A Cckm « (e/,A + DA^x - £k,A) /km C km- (10) 

Here, we introduced the local / energy, 

E fl x = Hf,x-DAx, (11) 

and defined D = l-(n{) + (n{)/vf and Ax = ^Z)k^k,A- 
The factors D in Eqs. {5J and (|10|) are caused by the Hub- 
bard operators in the renormalization ansatz Q8J. Simi- 
lar expressions without factors D have also been found 
in Ref. |10j where a renormalization ansatz consisting of 
fermionic quasi-particles has been used. 

As one can see from Eq. the operator product 

/km c k™ can also be interpreted as an approximate eigen- 
operator of the Liouville operator Lq x- The correspond- 
ing eigenvalues are excitation energies and can be used to 
rewrite 7Yi,a, 

Wl,A = Y &k >* Fk (/km C km + h ' C -) . 
k,m 

where the O functions 

k ,A = O (A - \s f ,x + DA^x - £k,A|) 

restrict the particle-hole excitations to transition energies 
smaller than A. 
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3.2 Renormalization equations 

Next we want to follow the discussion of Ref. ^U] to de- 
rive renormalization equations for the parameters of the 
renormalized Hamiltonian 7i\, It turns out that the ac- 
tual calculations are only slightly modified by the new 
renormalization ansatz JHJ which now includes correlated 
Hubbard operators. 

To evaluate the new Hamiltonian Tirx—AX) according 
to Eq. an unitary transformation has to be performed 
to eliminate excitations within the energy shell between 
(A — AX) and A. As in Ref. ^Ojj we use the following 
operator ansatz for the generator Xx,a\ of the unitary 
transformation, 



X 



k (A,ZiA) A k (X,AX) (fl m c kr , 



C km/ k ' 







where the <9 k (A, AX) are products of two functions, 



e k (x,AX) = e KX [i-o 



k,(A-/lA)J 



Note that the <9 k (A, Z\A) confine the excitations which 
have to be eliminated by the renormalization step from 
A to (A — Z\A). The unknown parameters A k (A, AX) have 
to be fixed in such a way so that only transition with en- 
ergies smaller than the new cut-off (A — AX) contribute to 

H(\-AX)- 

As described in Ref. ^H], equations for the parameters 
A k (A, AX) of the generator of the unitary transformation 
as well as for the parameters of the renormalized Hamil- 
tonian TL\ can be found by comparing the coefficients of 
the operators in the renormalization ansatz (jSJ at cutoff 
(A — AX) and in the explicitly evaluated unitary transfor- 
mation (JSJ). 

Thus, we obtain the following equations: 



A k (X,AX) = 

= G k (X,AX) 
2VD 



(12) 



arctan 



2VDV k 



Hf,x + D(A ktX -Ax) -£k,A 



£k,(A-zlA) — £k,A — 



1 



(13) 



- [/i/,A + D (A K x - Ax) - £k,A] 
x jcos 2VDA k (X,AX) -lj 
v'WkSin \2VDA k (X,AX)} , 



M/,(A-4A) - M/,A = 
= ~ n A? 53 [ £ k,(A-zlA) - £k,\] 



(15) 



k 

X 



V S ~ !) ( C k m C k 



Vf-1 1 

4£>3/2 N 



^ { [M/,A + D (Ac, A - A) - £k,A] 



sin j2-\/DA k (A, A\) 



2V.DVk ^ cos 



2\/ZL4 k (A,Z\A)] - l}} 



x ( /L c km + h.c. 



"/-I I 
2L> JV 



53 [m/,a - D (^k,A - ^a) - e k , A ] 



x A k (A, Z\A) (/^ m Ck m + h.c. 



(A-4A) 



#A 



(16) 



[M/,(A-4A) - M/,a] 

53 [ £ k,(A-ZiA) — £k,A] 



Note that besides the factor 1/D in Eq. I|14|l these renor- 
malization equations exactly agree with those derived in 
Ref. 10 . However, the underlying Hamiltonians differ sig- 
nificantly because now the renormalization ansatz JSJl con- 
tains correlation effects by means of the Hubbard opera- 
tors. It will turn out that the Hubbard operators not only 
complicate the further evaluation of the renormalization 
equations but also successfully prevent the system from 
unphysical multiple occupation of the / sites. 

In deriving the renormalization equations i|12|) - i|ltj|) a 
factorization approximation has been employed so that 
the obtained equations still depend on expectation values 
which have to be determined simultaneously (see Ref. [TU] 
for details). Furthermore, an expansion in 1/vf has been 
avoided (and spin fluctuations have been neglected) so 
that the derived renormalization equations are valid for 
large as well as small degeneracies Vf. The limit A — * 
provides the parameters e k , p,f, Z\ k , and E of the effective 
Hamiltonian H = Hx^o = Ho.x^o, 



n = 53 £k ci m c km + A/ 53 /km/i 

k,m k , m. 



km 



(17) 



k.m 



A 



k,(A--4A) 



ik.A 



1 

D 



we are interested in. Here, it is important to notice that 
the renormalized Hamiltonian 7i no longer contains the 
hybridization between conduction and localized electrons. 
However, Ti. is not a non-interacting fermionic system be- 
cause H still takes into account electronic correlations by 
means of the Hubbard operators . Note that these cor- 
L £ k,(A-ziA) — £ k,Aj (,14) relations turn out to be crucial for a description of integral 
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valent states. On the other hand, the Hubbard operators 
/km a ^ so cause challenging difficulties in the further the- 
oretical treatment because they do not obey the usual 
fermionic anticommutator relations. 



3.3 Approximate solutions 

In the following we want to develop a strategy to solve the 
renormalization equations i|12|) - i|16[) approximately. Here, 
similar approximations as in Ref. jlOj shall be used to de- 
couple the renormalization of the different k values. In 
this way, all relevant quantities can be expressed as func- 
tions of a renormalizcd / energy if which is determined 
by numerical minimization of the free energy. 

As in Ref. jlOj , we use the following approximations for 
further evaluation of the renormalization equations Q12jl- 
fTF 



(i) All expectation values (which occur due to the ex- 
ploited factorization approximation) are assumed to 
be independent from the renormalization parameter A 
and are calculated using the full Hamiltonian TL. 

(ii) To decouple the renormalization of the different k val- 
ues, the A dependence of the renormalizcd / level is 
neglected, [if.x — DA\ ~ £/• The spirit of this approx- 
imation is similar to that assumed in the SB theory 
where a renormalized / energy is also used from the 
very beginning. Note that if has to be interpreted as 
local / energy of the renormalized model (|17(l . 

At this point it is important to notice, that our old 
analytical solution of Ref. ^H] can be easily obtained if 
the Hubbard operators in the final Hamiltonian l|17(l are 
replaced by usual fermionic operators. Formally, one em- 
ploys 

( m ) Ekjn /km/ km ~ Sk,m /km/ km anC ^ 
(/km/km)NL ~ D (/ km /km)NL ■ 

to ensure that, on a mean-field level, the renormalized 
Hamiltonian does not generate unphysical states. How- 
ever, as already discussed above, the obtained effective 
model does not prevent anymore from multiple occupa- 
tion of / sites if (iii) has been employed. We have already 
argued that such an approximation can only lead to useful 
results as long as only very view / type states below the 
Fermi level are occupied. Thus, only HF-like solutions can 
be observed in this way. To obtain the analytical solution 
of Ref. JUj , one also has to employ 

for further simplification. 

In the following we only want to employ approxima- 
tions (i) and (ii). In particular, we keep the Hubbard op- 
erators in the final Hamiltonian 117|) so that both mixed 
valent and integral valent states can be described. 



Eqs. (|14(1 and (|16fl can be easily integrated between the 
lower cutoff A — > and the cutoff of the original model A, 



A v 



D 



[£k — £kj , 



(18) 



E = -N(h{) [if-e f ] + 



(19) 



As already mentioned above, the approximations (i), (ii) 
decouple the different k values from each other so that 
Eq. i|12[l and are completely similar to those obtained 
for the Fano- Anderson model (compare Ref. JOI)- Thus, 
two quasi-particle branches are obtained, 



- if +£k sgn(e / - e k ) 

£k = r VV k , 



(20) 



^ := i f + DA, = + + Sgn( V £k V k , (21) 



where 



Wv 



(e k -e / ) 2 +4^|y k |2. 



Note that the one-particle energies (I2U|I and 121|) still de- 
pend on two unknown quantities: the renormalized / level 
if and the / occupation number (h{) (that determines D 
as defined above). 

In Ref. ^U] all expectation values as well as the renor- 
malized / level if have been determined by functional 
derivative of the free energy. However, here, this approach 
can not easily be applied because the Hubbard operators 
contained in the renormalized Hamiltonian (|17J) do not 
fulfill the usual fermionic anti-commutator relations. Fur- 
thermore, the derivation of the free energy would also lead 
to problematic 6 functions that are caused by the abrupt 
change of the statistic of the quasi-particle excitations at 
if. (In Ref. ^H] these contributions do not appear because 
both c-like and /-like excitations are caused by fermionic 
quasi-particles.) Therefore, a different approach has to be 
developed to determine the renormalized / level if and 
the expectation values. 

In the following, the expectation values of the original 
Hamiltonian TL will be calculated using the renormalized 
one-particle operators as derived in Ref. |10| . 



/km 

where we defined 



°) =^kC k m + W km , 
0) = -Dv k ct m + u k fl 



km' 



|Mkl = 2\ 1 ~-^r sgn(e/ ~ ek) 

l^kl 2 = ttf:\ 1 + £k TTr £/ sgn(e~/ - e k ) 



2D 



(22) 
(23) 



(24) 
(25) 



Thus, the required expectation values of the full Hamil- 
tonian H can be traced back to those calculated with re- 
spect to the renormalized Hamiltonian TL because (A) = 
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lim A ^ (^(A))w A holds, 



C km C kn 



1 TT^-sgn (e/ - e k ) 



1 + 



£k — if 



/(£k) (26) 



sgn (if - e k ) 



/(«k) 



/km c km 



h.C. 



(27) 



= - 2 sgn (e, - £k ) M [ /(fk) _ / (Ck ) 
Here, we introduced the Fermi function 



/(£k) 



km ckm/ w - 1 + e/a g k ■ 



and defined 



/(Wk) 



n Ukr. 



./km) 



(28) 



Note that the factor D in Eq. (|28|) has been introduced to 
underline the similarities of Eq. (|26|l and 1(27(1 with the cor- 
responding results of the analytical treatment of Ref. |10| . 

In principle, the / occupation number (h{) could also 
be calculated using the renormalized one-particle opera- 
tors. However, here we alternatively employ the particle 
conservation under unitary transformations. Thus, we ob- 
tain 



(h{) 



(29) 



\Vf 



2 N 



E 



1+ — • / sgn(e / -£ k ) 



2 AT ^ 



/(Wk)- 



For actual calculations one needs to evaluate Eq. 1(28(1 
in order to determine the expectation values of the full 
Hamiltonian as given in Eqs. HU), (EH, and Because 
of the unusual properties of the Hubbard operators, there 
is no straightforward way to evaluate Eq. 128fl and further 
approximations are necessary. As long as the renormalized 
/ level is situated above the chemical potential a mean- 
field treatment of the electronic correlations contained in 
TL might be sufficient, and we would find /((2> k ) « f(&k) as 
directly obtained by employing approximation (iii) men- 
tioned above. On the other hand, here we are also inter- 
ested in solutions of the PAM with a renormalized / level 
below the Fermi level which require a theoretical treatment 
of the electronic correlations in TL beyond a mean-field ap- 
proximation. Therefore, Eq. 1)28(1 is evaluated as follows 



(/km /km) -fi 



/(^k) ( { /km ' /km j 



Tr e 



km e 



where the approximated / excitation energy as derived in 
Eq. I |l Up has been used. Thus, (|28|l can be rewritten as 



r/(wk) 



l 



■ Ek' /(^k' 



(30) 



Unfortunately, approximation ((30(1 does not offer a direct 
link to the mean-field result, /(cl>k) ~ /(wk), for renormal- 
ized / energies above the Fermi level. Thus, differences be- 
tween the presented treatment and the analytical solution 
of Ref. JU| will appear. 

At this point all physical quantities can be calculated 
as function of the renormalized / energy if. Because we 
have employed approximation (ii) it is not possible any- 
more to use the renormalization equation 1(15(1 for /j,f t \ to 
determine if. Therefore, the local / energy if is consid- 
ered as a free parameter and is determined by minimiza- 
tion of the free energy. Because of the unusual anticom- 
mutator relations of the Hubbard operators / krra , the free 
energy can not be directly determined. Instead, 



dF _ de k / t \ d/3/ / - t , 



di f ^ di f \ km "7h di f 



km 



k,m 



H 



^ dif ((^Mnl)^ + dif 

k,m 



(31) 



is numerically integrated in order to calculate the free en- 
ergy F as function of the renormalized / energy if. Note 
that Eq. I|31l) has been obtained from the renormalized 
Hamiltonian l|17|) . Actual results are discussed in the next 
section. 



4 Results 

It is believed that the one-particle energy Ef of the local- 
ized / electrons is smoothly changed in CeCu2Si2 due to 
pressure 0]. Therefore, we want to discuss the physical 
properties of the PAM as a function of ef. 

At first let us consider an one-dimensional PAM with a 
linear dispersion relation e k for the conduction electrons in 
the energy range between —1 and 1, and a k independent 
hybridization y k = V. The other parameters are chosen 
as follows VfV 2 = 0.36, chemical potential /i = 0, and 
T = 0.00001 where all energies are given in units of the 
half bandwidth. 

As one can see from Figs.n an d|3 we obtain two dif- 
ferent types of solutions depending on the value of the un- 
renormalized / level £/. First of all, we obtain the usual 
SB type solutions with intermediate valence states rif < 1 
where the renormalized energy if is energetically located 
above the Fermi energy. If the unrenormalized energy Ef 
is lowered the renormalization contributions are no longer 
sufficient to push if above the Fermi level, and the renor- 
malized / energy if is located far below the Fermi energy. 



A. Hiibsch and K. W. Becker: Valence transition in the periodic Anderson model 



7 



-0.8 




0.4 



» » this work 

Ref. 10 

■ ■ ■ ■ SB mean-field 



-0.8 -0.6 -0.4 -0.2 

7 



-0.4 



-0.8 



0.8 



0.6 



0.4 



~^ I 1 I 
(a) v,= 20 



+ 



— • this work 

■ - Ref. 10 

■ ■ ■ SB mean-field 



-0.8 




-0.6 



7 



Fig. 1. Renormalized / level if [panel (a)] and averaged / 
occupation number n/ = (n{) [panel (b)] as function of the 
unrenormalized / energy Sf where an one-dimensional PAM 
(N = 10000, v s = 4, v;V 2 = 0.36, fi = 0, T = 0.00001) with a 
linear dispersion relation for the conduction band in the energy 
range between —1 and 1 has been considered. (All energies 
are given in units of the half bandwidth.) For comparison, the 
results of the PRM approach of Ref. |10[ and of the SB mean- 
field theory are drawn with dashed and dotted lines. 



Fig. 2. Renormalized / level if [panel (a)] and averaged / 
occupation number n/ = (n{) [panel (b)] as function of the 
unrenormalized / energy ef for an one- dimensional PAM with 
Vf = 20. Other parameters are chosen as in Fig. 



because both do not explicitly ensure that / sites can only 
be either empty or singly occupied as already discussed 
above. 



In this case, the averaged / occupation nf is almost ex- 
actly 1 and an integral valence state is obtained. 

Figs. Hand also reveal the very good agreement be- 
tween the HF type solutions of the presented PRM ap- 
proach and the analytical results of Ref. 10j ■ In this way 
it is proven that the Hubbard operators can be replaced 
by usual fermionic operators (compare approximation (iii) 
in Sec. l3.3fl because in this case only very few / type states 
below the Fermi level are occupied as discussed above. In 
this regard one needs to keep in mind an important differ- 
ence between the SB theory and our PRM approach: The 
quasi-particles of the SB theory change their character as 
function of k between more /-like and more c-like behav- 
ior. In the PRM excitations do not change their charac- 
ter as function of k, and the quasi-particle energies show 
jumps in their k dependence if £/ is energetically located 
within the conduction band. Note, however, that the vari- 
ous parts of the quasiparticle bands fit perfectly together, 
as one can see from Eqs. (l%0l and ifSTjl. 

For comparison, the results of the analytical solution of 
Ref. lij and of the SB theory are shown as well in Figs.H 
and [21 As one can see, no solution with renormalized / 
level if < could be found for these analytical approaches 



The well-defined transition between the two different 
solution types is of particular interest. As expected, for 
the HF-like solution the /-charge is always smaller than 1 
due to hybridization processes between / and c-electrons. 
Simultaneously heavy quasiparticle bands are formed at 
the Fermi surface. To describe the HF behavior the full 
Anderson model has to be considered. As the bare /-level 
moves to smaller energies a transition to an integral va- 
lence charge of rif = 1 is observed (similar to the Ander- 
son impurity model |14j). In this case only the c electrons 
should form the Fermi surface. Thus, the observed valence 
transition can also be interpreted as a collapse of the large 
Fermi surface of the HF state which is formed by conduc- 
tion as well as by localized / electrons. Note, however, that 
the question whether localized electrons contribute to the 
Fermi sea volume or not is still controversially discussed 
in the literature |15| . 

As one can see from Figs.^andEJ the obtained valence 
transition is much more pronounced for small degenera- 
cies Vf, and a smooth transition can be expected in the 
limit Vf — > oo of the SB theory. Therefore, a sharp valence 
change in generalized SB theories can only be obtained if a 
rather large additional Coulomb repulsion between / and 
conduction electrons is present in the system How- 
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Fig. 3. Panel (a) and (b) [(c) and (d)] show the results 
for the renormalized / level if and the / occupation number 
rif = (h{) for a two-dimensional [three-dimensional] system 
with 100 x 100 [26 x 26 x 26] lattice sites. As in Figs. [T]and|l 
a linear dispersion relation £k = e(|k|) has been chosen, and 
the electronic band covers an energy range between -1.5 and 
0.5 [-1.8 and 0.2] where n = 0, v f = 4, u f V 2 = 0.36, and 
T = 0.00001. (Energies are given in units of the half band- 
width.) Note that the conduction band has been energetically 
shifted in in order to ensure a filling comparable to the one- 
dimensional case of Fig. because a smaller filling of the con- 
duction band reduces the change in the / occupation n / at the 
transition point. 

ever, here we have shown that such a valence transition 
can also be obtained in the plain PAM if corrections for 
small degeneracies Vf are properly taken into account. 

One of the advantages of the analytical PRM is the 
opportunity to consider much larger systems than acces- 
sible by numerical methods. Therefore, we are also able to 
study two- and three-dimensional systems of reasonable 
sizes. In this way we can easily show that the observed 
valence transition is not an unique phenomenon of the 
one-dimensional PAM. The valence transition also occurs 
in two- and three-dimensional systems as can been seen 
in Fig. |3 Therefore, the observed behavior has to be con- 
sidered as a general feature of the PAM, and our results 
should also be of relevance for actual physical HF systems 
like CeCu2Si2 or related compounds. 



5 Discussion and Summary 

The occurrence of a valence transition in the plain PAM is 
the main finding of this paper. In contrast, a rather large 
additional Coulomb repulsion has been claimed to be nec- 
essary for the valence transition in an extended PAM H 



IHj. The studies of Refs. 0]|H] were based on a slave-boson 
fluctuation approximation that extends the well-known 
slave-boson mean-field theory |11II12| but still employs the 
limit of large degeneracy Vf — > oo. Our results show (com- 
pare Figs. and [2| that the observed valence transition 
becomes smooth in this limit. Therefore, it is reasonable 
that an additional interaction was found to be necessary 
in order to obtain a valence transition in an approach em- 
ploying Vf — > oo. 

Our work also shows the importance of taking care of 
a physical / occupation in theoretical approaches. In par- 
ticular, it turns out that a completely uncorrelated model 
is not able to prevent from unphysical multiple occupa- 
tion of / sites, and no integral valence states can be found 
in this way. In contrast, the presented PRM approach to 
the PAM explicitly suppresses unphysical multiple / occu- 
pation which is, in particular, crucial for integral valence 
states. 

We obtain two solution types: a mixed valence state 
with a renormalized / level s / above the Fermi energy and 
an integral valence state with if below the Fermi level. 
Furthermore, parameter regimes exist where the transi- 
tion between the two solution types is accompanied by a 
drastic change in the / occupation. Such a sharp valence 
transition occurs in one-dimensional as well as in two- and 
three-dimensional systems so that this behavior has to be 
considered as a general feature of the PAM. Note that a 
similar valence transition has been experimentally found 
in CeCu2Si2 from high pressure experiments 

In the case of an integral valence state one would ex- 
pect that the system can be described by a Kondo Hamil- 
tonian which is gained from the PAM by the Schrieffer- 
Wolff transformation 16 for F/|e/| <C 1. Note, however, 
that in the present approach spin fluctuations have been 
neglected altogether but Kondo-like and RKKY-like in- 
teractions as well as higher charge fluctuation terms are 
automatically generated during the renormalization pro- 
cedure. These contributions will have to be considered in 
the future. One might expect that additional spin and 
charge fluctuations might possibly give rise to magnetic 
and superconducting phases both for the intermediate va- 
lence and for the integer valence regime. Also, one may 
speculate that the magnitude of the magnetic moment 
will be different for these cases due to additional screen- 
ing processes. 

The PRM approach presented in this paper only ad- 
dresses the question for the valence transition in the plain 
PAM. However, as mentioned above, the PRM scheme of- 
fers great opportunities to include additional interactions 
which are automatically generated during the renormal- 
ization procedure. Therefore, extensions of the PRM treat- 
ment might be promising starting points to study the com- 
peting interactions in CeCu2Si2 and related compounds in 
more detail. 
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